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ABSTRACT 

We use a nearly complete sample of Gamma Ray Bursts (GRBs) detected by the Swift satellite 
to study the correlations between the spectral peak energy iSpeak of the prompt emission, the 
isotropic energetics Ei^o and the isotropic luminosity Liso- This GRB sample is characterized 
by a high level of completeness in redshift (90%). This allows us to probe in an unbiased way 
the issue related to the physical origin of these correlations against selection effects. We find 
that one burst, GRB 061021, is an outlier to the fipeak - Emo correlation. Despite this case, we 
find strong fipeak - ^iso and /Speak - iiso correlations for the bursts of the complete sample. 
Their slopes, normalizations and dispersions are consistent with those found with the whole 
sample of bursts with measured redshift and /ipeak- This means that the biases present in the 
total sample commonly used to study these correlations do not affect their properties. Finally, 
we also find no evolution with redshift of the E^^nk - Ei^o and /Speak - i-iso correlations. 
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1 INTRODUCTION 

, One of the most debated issues in Gamma Ray Bursts (GRBs) con- 
'= cerns the existence of correlations among the spectral parameters of 
the prompt emission and its energetic and luminosity. Three robust 
correlations have been identified. Each of them involves the rest 
frame peak energy iipeai^ of the vFy prompt spectrum. This quan- 
tity strongly correlates with (i) the isotropic energy E^^a (Amati re- 
lation, Amati et al. 2002), (ii) the isotropic peak luminosity Liso 
(Yonetoku relation, Yonetoku et al. 2004) and (iii) the collimation- 
corrected energy Ey (Ghirlanda relation, Ghirlanda et al. 2004b; see 
also Liang & Zhang 2005). These correlations are valid for long 
GRBs. There are indications that short bursts (i.e. bursts with ob- 
served duration <2 s) obey the very same £peak - ^iso relation de- 
fined by long events but they are inconsistent with the iipeak - £^iso 
correlation (Ghirlanda et al. 2009; Amati 2010). 

The spectral-energy correlations have relevant implications 
both for the theoretical understanding of the burst physics and for 
the application of GRBs as cosmological tools (Ghirlanda et al. 
2004a; Ghirlanda et al. 2006). While there is no consensus yet 
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on the interpretation of these correlations (e.g. Eichler & Levin- 
son 2004; Levinson & Eichler 2005; Rees & Meszaros 2005; 
Thompson, Meszaros & Rees 2007), an interesting progress has 
been recently made by studying the comoving properties of GRBs 
(Ghirlanda et al. 2012; Panaitescu et al. 2009). 

The £peak - £^iso and iipeak - ^iso correlations have been dis- 
covered with a very small sample of 12-16 GRBs detected by 
BeppoSAX and/or BATSE (Amati et al. 2002; Yonetoku et al. 
2004). They have been subsequently confirmed by adding bursts 
detected by other satellites (e.g. HETE 11, KonusPNlND, INTE- 
GRAL, Swift and Fermi). At present the sample of bursts defining 
the Speak - £iso ^ud Speak - ^iso Correlations comprises more than 
130 events, with only two well known (very peculiar) outliers (but 
see Ghisellini et al. 2006). 

Despite the considerably large number of bursts consistent 
with these correlations, their physical origin is debated. Some au- 
thors claim that they are the result of instrumental selection ef- 
fects (Nakar & Piran 2005; Band & Preece 2005; Butler et al. 
2007; Shahmoradi & Nemiroff 2011). Other studies (Ghirlanda 
et al. 2008; Nava et al. 2008) quantified the possible instrumen- 
tal selection biases finding that, even if they do affect the sam- 
ple, they cannot be responsible for the existence of the spectral- 
energy correlations. Ghirlanda et al. (2010; 2011a; 2011b) added 
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an important piece of information to this debate. They analysed the 
time-resolved spectra of both long and short bursts detected by the 
Fermi satellite. They found that the time-resolved E^.^ correlates 
with the luminosity during the temporal evolution of the burst (and 
for bursts without z they derive a correlation between the time- 
resolved E'^^j^ and the flux). This is a feature common to short and 
long GRBs. The time-resolved correlations found within individual 
GRBs are similar to each other and they are similar to the £'peak-^'iso 
relation defined by the time-integrated spectral properties of differ- 
ent bursts (i.e. the Yonetoku relation). This is a strong argument in 
favour of the physical origin of the correlations because it is very 
unlikely that within a single burst the instrumental selection effects 
(e.g. the detector threshold) or possible effects related to the depen- 
dence of the peak energy, energetic and luminosity on the redshift 
can play any role. 

GRBs added to the £pcak-£^iso and -Epeak-^iso correlations need 
to have their redshifts measured. This raised the suspect that these 
correlations might be biased by the high level of incompleteness in 
redshift of the samples of GRBs so far used. In this paper we study, 
for the first time, the Epeak - £^iso and Epc^k - ^^iso correlations with 
a sample of GRBs with a high level of completeness (90%). The 
sample is presented by Salvaterra et al. (2012) and it comprises 58 
events detected by Swift with a 1-s peak photon flux P > 2.6 ph 
s"'cm"^ (integrated in the 15-150 keV energy range). Fifty-two 
out of 58 GRBs have measured redshift z, while in 3 cases an upper 
limit on z can be set from the photometry of the afterglow or of the 
host galaxy. The completeness level then increases up to ~95% by 
considering the redshift constraints imposed by the detection of the 
afterglow or host galaxy in some optical filters. 

The complete sample used in this work is particularly ap- 
propriate to study the possible evolution of the E^^.j^ - E-^so and 
Epeak - E\so correlations with redshift. Li (2007), considered a sam- 
ple of 48 GRBs (from Amati 2006, 2007) and investigated if the 
£peak - Emo cotTelatiou evolves with redshift, finding that it does, 
becoming steeper at higher redshifts. Ghirlanda et al. (2008) re- 
peated the same test with a larger sample of 76 bursts and found no 
evidence of evolution with redshift. These results, however, were 
obtained by considering highly incomplete GRB samples. 

Our goals are (i) to study how Speak is correlated to E\sa and 
Liso within the complete sample of bursts and (ii) to study the pos- 
sible evolution of these correlations with redshift through this unbi- 
ased sample. In §2 we present the spectral properties of the GRBs 
composing the complete sample. We study the £peak - £^iso and 
£peak - iiso Correlations defined by this sample in §3 and compare 
their properties with those defined by the (incomplete) larger sam- 
ple of GRBs with measured z and E'^^. In §4 we discuss and sum- 
marise our conclusions. Throughout the paper we assume a stan- 
dard cosmology with h = = 0.7 and fl,,, = 0.3. 



2 THE SAMPLE 

In order to study a relatively unbiased redshift distribution of 
GRBs, Jakobsson et al. (2006) carefully selected a sample contain- 
ing bursts which have 'observing conditions' favourable for red- 
shift determination (see their web pag^for an updated compilation 
of the sample and an explanation of the selection criteria). This rel- 
atively 'clean' sample has a redshift recovery rate of roughly 50%. 
In order to reach a higher level of completeness, Salvaterra et al. 
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(2012) considered this sample and also required that the photon 
peak flux, measured by Swi/r/BAT in the 15-150 keV energy band, 
is f > 2.6 ph s^'cm"-. 

Up to May 2011, 58 GRBs match this selection criteria and 
are listed in Table [T] Fifty-two of them have measured redshift so 
that our completeness level is 90%. Of these 52, all but two (namely 
GRB 070521 and GRB 080602) have spectroscopic confirmed red- 
shift either from absorption lines over-imposed on the GRB optical 
afterglow or from emission lines of the GRB host galaxy. Moreover, 
for 3 of the 6 bursts lacking measured z the afterglow or the host 
galaxy have been detected in at least one optical filter, so that ~95% 
of the bursts in this sample have a constrained redshift. While this 
sample represents only ~10% of the full Swift sample, it contains 
more than 30% of long GRBs with known redshift. 

The bursts belonging to this sample (updated to May 2011) 
are listed in Table[T] We will refer to this sample as the 'complete' 
one. 

For all bursts in this sample, we collected from the literature 
the prompt emission spectral parameters and computed the rest 
frame peak energy Epc^=E°'^^^^(l + z), the isotropic energy Ei^o and 
the isotropic peak luminosity Lis^. The two latter quantities have 
been calculated in the rest frame energy range 1 keV-10 MeV. The 
narrow sensitivity range of the BAT instrument allows to determine 
^peak °'^ly when it falls approximately between ~ 15 keV and ~150 
keV. This condition is necessary but not sufficient: as demonstrated 
in the spectral simulation studies by Sakamoto et al. (2009), even 
when E'^^^, lies in the BAT energy range, the value of E°'^^^_ can be 
determined only if the spectrum is characterized by a high signal- 
to-noise. Similar conclusions were reached through numerical sim- 
ulations by Ghirlanda et al. (2008) and Nava et al. (2008). Almost 
all GRBs of the complete sample, however, have been detected 
also by other instruments with wider sensitivity ranges (mainly 
KonusfWWT) and, for bursts after August 2008, Fermi). This al- 
lowed to estimate E°'^'^^_ for a large fraction of events. We, there- 
fore, obtained a sample of 46 bursts (out of 58) with £peak, ^iso, 
Liso and redshift z, which are the quantities required to investigate 
the £peak - E\so and Speak - Siso correlations. For 6 GRBs there is 
a good knowledge of the prompt spectral shape (i.e. Speak is well 
constrained), but the redshift is unknown (in 3 cases it is possible 
to derive an upper limit on z from the photometry of the afterglow 
or of the host galaxy). For the remaining 6 bursts the determination 
of SpJ^,^ is not possible: in 3 cases the spectrum is well described by 
a single power-law, while in the other 3 cases it is possible to set 
a lower limit to the observed value of S"^^,,. Redshifts and spectral 
properties of the 58 GRBs of the complete sample are reported in 
Table [T] Here we also report the spectral parameters of the prompt 
emission spectrum (with their uncertainties) and the fluence and 
peak flux (with the energy ranges over which they are computed). 
These are the informations used to compute the rest frame isotropic 
energy and luminosities of the bursts in our samples. We also report 
for each burst the instrument used to derive the spectral parameters 
listed in the table. To summarize, our sample is composed as fol- 
lows: 

• 46 GRBs have measured z and Speak (and consequently E^so 
and Liso); 

• 3 GRBs have measured Speak but only an upper limit on z; 

• 3 GRBs have measured Speak but no redshift estimate; 

• 3 GRBs have measured z but only a lower limit on E°^^^_; 

• 3 GRBs have measured z but no S°^^ estimate; 

It is interesting to compare the Speak - Ei^o and Speak - Sjso cor- 
relations defined by the complete sample with the same correlations 
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defined with a larger sample of GRBs, comprising all bursts with 
measured redshifts and E^^^^ detected by different satellites. We up- 
dated also this sample (called in the following the 'total' sample) to 
May 2011. The total sample contains 136 bursts, for one of which 
the peak luminosity cannot be derived, due to the lack of measured 
peak flux. For GRBs detected by the Gamma ray Burst Monitor 
(GBM - Meegan et al. 2009) onboard Fermi after March 2010 we 
retrieved the public data from the online catalog and we analysed 
their time-integrated spectra following the standard procedure de- 
scribed in Nava et al. (2011a). This sample is based on the one 
presented in Ghirlanda et al. (2012) to which we added four bursts 
(GRB 060306, GRB 061021, GRB 081007, GRB 081221 belong- 
ing also to the complete sample - see TableO. 

We can thus compare the spectral-energy correlations defined 
by the complete and by the total samples. To point out possible 
differences we also consider the sample obtained from the total one 
by excluding the bursts of the complete sample (we will refer to it 
as the 'complementary' sample). 



3 RESULTS 

We first test the E^^.j^ - E\^a and the £peak - ^iso correlations with 
the sample of 46 bursts (included in the complete sample) with firm 
estimates of the redshift and of the spectral properties. The results 
are shown in Fig. [T] Data points from the total sample are plot- 
ted for comparison as grey dots. Red squares show the bursts of 
the complete sample. We estimate the Spearman's rank correlation 
coefficient and the associated chance probability for the complete, 
total and complementary sample and report these numbers in Table 
[2l We model the distribution of the data points in the logarithmic 
plane of Fig.[T]with a linear function. We applied the ordinary least 
squares bisector method (Isobe et al. 1990) for fitting the data. This 
choice is motivated by the fact that a priori there is no reason to 
assume either Speak or £iso (or Liso) as the independent variable, 
and also by the large scatter of the data points. The errors on the 
slope and normalization are determined by fitting in the barycenter 
of data points, where they are uncorrelated. The results are given 
in Table[2l In Fig.[T]the best fit lines are shown for each sample as 
follows: a solid black line for the total sample, a dashed red line for 
the complete sample and a dot-dashed black line for the comple- 
mentary sample. The shaded region marks the 3cr of the scatter dis- 
tribution of the data points around the best fit line for the total sam- 
ple. The scatter estimated for the other two samples is also listed 
in Table [2] In Fig. [T] we also report the position of GRB 980425 
and GRB 031203 (blue filled stars), which are well known outliers 
of both the Epc^k - £iso and the Epc^k - ^^iso correlation. These two 
bursts have peculiar properties and they are generally considered as 
outliers to the confirmed relations. 

The behaviour of GRB 061021 is very peculiar, since it lies at 
more than the 3cr from the Amati and Yonetoku correlations (con- 
sidering the scatter defined by the total sample - shaded regions in 
Fig.B- 

Twelve events belonging to the complete sample cannot be di- 
rectly used to test the correlations, since they have some unknown 
property: the redshift or the spectral peak energy (see Sec. 2). How- 
ever, we can still test their consistency with the correlations (results 
are shown in Fig.[2j. To check the consistency different methods are 
adopted for the following different cases: 

• lower limit on E'^^^ and secure estimate of z (3 cases): we 
derive the lower limit on the rest frame peak energy and on E-^^^ and 
Liso (right and up arrows in Fig.|2li; 



• no estimate of E'^^^ and secure estimate of z (3 cases): we 
consider the energy range of sensitivity of the instrument and the 
photon index of the power-law fit. Photon indices greater (lower) 
than -2 [in the notation N{E) oz E'^] suggest peak energies greater 
(lower) than the upper (bottom) energy edge of the instrument. An 
upper or lower limit to Speak (and consequently on Ei^o and Lj^o) can 
be set (empty triangles in Fig.|2l; 

• spectral properties well constrained but upper limit on z (3 
cases): Speaks ^iso and Ljso are estimated for a variable redshift rang- 
ing from 0.01 to its upper limit (coloured continuous curves in Fig. 

IS; 

• spectral properties well constrained but no redshift estimate (3 
cases): Speak, £iso and Liso are estimated for values of the redshift 
ranging from 0.01 to 10 (coloured continuous curves in Fig.[2]l. 

3.1 Evolution with redsliift 

The possibility that evolutionary effects play a role in the spectral- 
energy correlations still represents an open issue. To investigate if 
the slopes of the Speak - £iso and Speak - hso correlations have a 
dependence on the redshift we repeat the same method used by Li 
(2007) and Ghirlanda et al. (2008) and apply it to our complete 
sample. We exclude from the fit the burst which lies at more than 
3cr of the Amati and Yonetoku relations. We define 4 sub-samples, 
with almost the same number of bursts, divided as follows: (i) z < 1 
(10 events, yellow empty circles in Fig. [3j, (ii) 1 < z < 1.7 (11 
events, filled purple stars), (iii) 1.7 < z < 2.4 (12 events, filled 
green triangles) and (iv) z > 2.4 (12 events, empty red squares). 
Bursts belonging to the total sample are not used for this test, since 
we want to take advantage of the completeness of our sample. The 
best fit slope of each sub-sample is reported in the inset of Fig. [3] 
both for the Amati (left panel) and Yonetoku (right panel) correla- 
tion. The light blue horizontal stripe, shown in the inset, represents 
the two-sided Icr interval of the slope derived by considering the 
whole complete sample and by excluding (for consistency) the two 
outliers. There is no evidence of a systematic evolution of the slope 
with z. The slopes defined by individual bins of redshift are all con- 
sistent, within Icr, with the slope defined by the entire complete 
sample. This result does not change by choosing a different num- 
ber of bins. 

Bursts with the highest redshifts populate only the high- 
energy (and high-luminosity) region of the Amati and Yonetoku 
correlations. This is due to the fact that we selected relatively bright 
bursts, by requiring a minimum peak flux. The effect of this selec- 
tion can be investigated as follows. We consider a series of spectra 

characterised by typical values a 1 and /3 2.3. We impose 

that the photon peak flux in the 15-150 keV energy range is equal 
to our limiting flux (P = 2.6 ph s^'cm"-). For different values of 
^pcik derive the corresponding limiting luminosity by consid- 
ering two cases: (i) z = 2 (dashed curve in Fig.O and (ii) z = 3 
(dot-dashed curve). We consider these redshifts since they are rep- 
resentative of the two sub-samples at higher redshifts, which popu- 
late only the high-energy (and high-luminosity) part of the planes. 
From Fig. [3] it is clear that this effect is due to the use of a flux 
limited sample. 

We may wonder if these limiting luminosities derived at dif- 
ferent redshifts may affect our study on the redshift evolution. We 
identify two possible effects: at high redshift the range of avail- 
able luminosities is smaller, i.e. the reduced leverage in luminos- 
ity might harden the estimate of the correlation slope and enlarge 
the associated uncertainty. Moreover, since the introduced cutoff is 
nearly vertical in the Speak-Siso plane, we expect a possible steep- 
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Eiso [erg] L^so [erg/s] 

Figure 1. Amati (left panel) and Yonetoku (right panel) correlations. Grey filled circles refer to the total sample. Their power-law fit is shown as a solid dark 
line. The shaded region represents the 3it scatter of the distribution of points around this best fit line. Bursts belonging to the complete sample (and with secure 
redshift and Speak estimate - 46 events) are marked as empty red squares and the red dashed line represents their best fit model. The dot-dashed line is the fit 
to the complementary sample. The position of the two historical outliers, GRB 980425 and GRB 031203, is also shown (blue stars), together with the burst of 
our complete sample, GRB 061021, that is above the 3it scatter limit. 



ening of the slopes toward high redshift. This effect has been inves- 
tigated also by Li (2007). By performing simulations, he concludes 
that a spurious evolution of the slope with the redshift can be in- 
troduced, with high-redshift bins defining steeper relations. We do 
not find signs of this effect in our data (see the inset in Fig.[3]l, indi- 
cating that, even if at high redshift we are missing low-luminosity 
events, we are still able to recover the slope of the correlations 
which are consistent with those derived with low-z bins. 



3.2 Distributions of £pcak) ^iso and L-^^o 

The complete sample adopted in this work has been well selected 
and the effects induced by the limiting peak flux have been well un- 
derstood and quantified in Sect. 3.1. On the other hand, up to now, 
the spectral-energy correlations have been studied by means of a 
sample which includes all possible GRBs with known redshift and 
Epeak, detected by different instruments. This total sample, there- 
fore, is affected by several instrumental effects that can play im- 
portant roles (as the different detector sensitivities and the different 
energy ranges on which the detection is performed). To estimate 
these effects and to understand how they influence the properties of 
the total sample is not trivial. Here we can use a different approach 
and point out possible biases by comparing the results obtained 
from the total sample with those from the complete one. From the 
point of view of correlations, we demonstrated that the total sample 
does not introduce evident biases on their slopes, normalizations 
and scatters (see Table |2ll. However, this does not necessarily im- 
ply that the complete and total samples share the same properties 
in terms of distributions of £peak, £iso and L\^o. 

We investigate possible effects present in the total sample by 
comparing the £pcak> £iso and Lisa distributions of the complete, to- 
tal and complementary sample. For each quantity, we also apply 
the Kolmogorov-Smirnov (K-S) test and estimate the probability 
that two distributions are drawn from the same parent population. 



In particular, we apply this test to compare i) the complete and to- 
tal sample and ii) the complete and complementary sample. Re- 
sults are shown in Fig.|4] In all cases, the K-S probability between 
the complete and the total sample is large, and the distributions are 
consistent within \a. The probability slightly decreases when com- 
paring the complete and the complementary sample. In this case all 
distributions are consistent within 2a. A slight effect is possibly 
present in the Ljso distribution, where the selection cut discussed in 
the previous section lowers the number of intermediate/low lumi- 
nosity events. Due to the relation between L-^^a and iipeak, this effect 
also influence the £peak distribution. Even if numbers are still small 
for the complete sample, our results suggest that the total sample is 
not affected by any relevant bias. 



4 DISCUSSION AND CONCLUSIONS 

The definition of a sample of GRBs with a high level of complete- 
ness (90%) allows to address several open issues regarding popula- 
tion studies of GRBs. In this paper we use this complete sample to 
study the spectral-energy correlations in an unbiased way. Out of 
58 GRBs satisfying the selection criteria 46 have well defined spec- 
tral properties and measured redshift. In particular, the rest frame 
peak energy £pcak can be determined and its relation with E\^a and 
Liso can be investigated. We find that this sample defines two strong 
correlations: (i) ZJpeak - Ei^o (rank's correlation coefficient p = 0.76 
and chance probability Pchance = 7 x 10"'") and (ii) iipcak - ^'iso 
(p = 0.65 and P,ha„ce = 1 x lO""). 

The slope of the correlations defined by the complete sam- 
ple are 0.61 and 0.53 for the Amati and the Yonetoku correlations, 
respectively. These slopes are consistent within Icr (2cr) with the 
slopes of the £peak - ^iso (£^peak " £^iso) Correlation defined by bursts 
not satisfying the completeness criteria (complementary sample). 
From Table |2] we also note that the scatter of the correlations de- 
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Ei,o [erg] Lj^o [erg/s] 

Figure 2. The consistency with the Amati (left panel) and Yonetoku (right panel) correlations of bursts with unknown redshift or E°^^ is shown. For bursts 
without redshift estimate (or with an upper hmit on the redshift) the test is performed by varying z from 0.01 to 10 (or to the upper limit). The resulting iJpeak. 
£iso and Ljso are shown as continuous lines. Different colours refers to different ranges of redshift (see legend). Arrows mark those bursts for which z is known 
but there is only a lower Umit on the E'^^^ value (which transforms into a lower limit on £iso and Liso). Triangles represent those events whose spectrum is 
well described by a single power-law model. In these cases we can set a lower limit (upwards triangles) on fipcak. ^iso and Liso by considering the photon 
index of the power-law fit and the energy range of sensitivity of the instrument. The total sample (and its 3(T scatter - shaded region) and the complete sample 
are also shown for reference. The position of the two histoiical outliers, GRB 980425 and GRB 031203, is also shown (blue stars). 
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Figure 3. Evolution with redshift for both the £peaic - Liso (left panel) and the fpeak ~ ^iso (right panel) correlation. GRBs from the complete sample have 
been divided into 4 sub-samples, accordingly to their redshift (see legend). The inset shows the slope of the correlation for each sub-sample. The shaded blue 
region is the two sided Icr interval of the slope found by considering the whole complete sample of 46 GRBs and by excluding for consistency GRB 061021. 
Big symbols mai'k the centre of mass of data points for each subsample and for the whole sample (dark cross). Grey filled circles show data points from the 
total sample. The dashed black curve marks the lowest possible luminosity of the flux limited sample estimated for GRBs with z = 2. The dot-dashed curve 
represents the same for bursts with z = 3. 
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Table 1. List of bursts that define the complete sample. Redshifts (from Salvaterra et al. 2012), spectral parameters, rest frame peak energies fipeak. isotropic 
energies Ej^ and luminosities Lj^o are provided. £iso and Lj^ are estimated in the l-l(f keV energy range. All quoted errors are at the 90% confidence 
level. Column 8 reports the name of the mission from which spectral properties have been derived: S=Swift, K=Konus/Wind, F=Fermi, Su=Suzaku, H=Hete. 
References in the last columns are for the spectral properties: [1] Perri et al. 2005, [2] Sakamoto et al. 2008, [3] Golenetskii et al. 2005a, [4] Sakamoto et al. 
2006, [5] Blustin et al. 2005, [6] Butler et al. 2007, [7] Sakamoto et al. 201 lb, [8] Crew et al. 2005, [9] Golenetskii et al. 2005b, [10] Cabrera et al. 2007, [11] 
Sakamoto et al. 2006b, [12] Barthelmy et al. 2006, [13] Amati et al. 2007b, [14] Golenetskii et al. 2006a, [15] Sakamoto et al. 2011, [16] Golenetskii et al. 
2006b, [17] Golenetskii et al. 2()06c, [18] Golenetskii et al. 2()07a, [19] Golenetskii et al. 20()7b, [20] Krimm et al. 2007, [21] Golenetskii et al. 20()7c, [22] 
Racusin et al. 2008, [23] Golenetskii et al. 2()08a, [24] Golenetskii et al. 2008b, [25] Golenetskii et al. 2008c, [26] Golenetskii et al. 2008d, [27] Golenetskii et 
al. 2008e, [28] Golenetskii et al. 2008f, [29] Starling et al. 2009, [30] Nava et al. 2011a, [31] Bissaldi et al. 2008, [32] Golenetskii et al. 2008g, [33] Ukwatta 
et al. 2008, [34] Golenetskii et al. 2009a, [35] Golenetskii et al. 2009b, [36] Golenetskii et al. 2009c, [37] Pal'shin et al. 2009, [38] Noda et al. 2009, [39] 
Golenetskii et al. 2009d, [40] Foley et al. 2010, [41] Golenetskii et al. 2010, [42] this work, [43] Golenetskii et al. 2011a, [44] Golenetskii et al. 2011b. 
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(2.2±0.17)xl0"^ 


15-150 


S 


115±27 


1.69±0.17 


4.76±0.37 


1, 2 


050401 


2.90 


-1.0 [-2.45] 


19.3±0.4 


20-2000 


(2.45±0.12)xl0"* 


20-2000 


K 


499±117 


40.6±0.84 


201±9.85 


3 


0504 16A 


0.653 


-1.0 [-3.4] 


0.35±O.O3 


15-150 


5.0±0.5'' 


15-150 


S 


26±4 


0.094±0.008 


0.97±0.09 


4 


050525A 


0.606 


-0.99±0.11 


20.1±0.5 


15-350 


47.7±1.2° 


15-350 


S 


127±6 


2.32±0.06 


7.23±0.18 


5 


050802 


1.71 


-1.6±0.1 


2.7±0.4 


15-350 


(22.1±3.53)xl0"* 


15-150 


S 


>192 


>3.19 


>8.90 


6,7 


050922C 


2.198 


-0.83±0.24 


3.14±0.31 


2-400 


(4.5±0.7)xl0"'' 


20-2000 


H/K 


416±118 


3.74±0.37 


184±28.7 


8,9 


060206 


4.048 


-1.12±0.30 


0.84±0.04 


15-150 


(2.02±0.13)xl0-'' 


15-150 


S 


409±116 


4.10±0.21 


49.6±3.24 


7 


060210 


3.91 


-1.12±0.26 


6.92±0.37 


15-150 


2.8±0.3° 


15-150 


s 


574±187 


35.3±1.9 


52.8±5.66 


10,11 


060306 


3.5 


-1.2 ±0.5 


2.5±0.34 


15-350 


(47.1±2.78)xl0"^ 


15-150 


s 


315±135 


8.26±1.12 


83.0±4.9 


6,7 


060614 


0.125 




21.7±0.4 


15-150 


11.6±0.7'' 


15-150 


s 


55±45 


0.25±0.10 


0.053±0.014 


12,13 


060814 


1.92 


-1.43±0.16 


26.9±2.6 


20-1000 


(2.13±0.35)xlO"* 


20-1000 


K 


750±245 


30.7±2.97) 


70.9±11.7 


14 


060904A 




-1.22±0.05 


8.49±0.50 


15-1000 


(1.3±0.3)xUr* 


20-10000 


K/Su/S 


235±25* 






15 


060908 


1.88 


-0.93±0.25 


2.81 ±0.11 


15-150 


(2.81±0.23)xl0"'' 


15-150 


S 


426±207 


4.41±0.18 


12.7±1.04 


7 


06091 2 A^ 


0.94 


-1.85±0.08 


2.87±0.42 


15-1000 


(2.5±0.9)x Ur'' 


20-10000 


K/Su/S 








15 


060927 


5.47 


-0.81±0.36 


1.12±0.07 


15-150 


(2.47±0.17)xl0"' 


15-150 


S 


459±90 


7.56±0.46 


108±7.6 


7 


061007 


1.261 


-0.75±0.02 [-2.79±0.09] 


193±2.62 


15-1000 


(1.2±0.1)xl0-^ 


20-10000 


K/Su/S 


965±27 


101±1.4 


109±9.1 


15 


061021 


0.346 


-1.22±0.13 


13.4±2.3 


20-2000 


(3.72±0.93)xl0"* 


20-2000 


K 


1046±485 


0.46±0.08 


1.73±0.43 


16 


061121 


1.314 


-1.32±0.05 


56.7±3.9 


20-5000 


(1.28±0.17)xl0-5 


20-5000 


K 


1402±185 


27.2±1.87 


142±18.9 


17 


061222A 


2.09 


-1.00±0.05 [-2.32±0.38] 


16.2±0.68 


15-1000 


(4.8±1.3)xl0"* 


20-10000 


K/Su/S 


1091±167 


22.5±0.94 


140±38 


15 


070306 


1.50 


^1.67±0.1 


9.0±0.63 


15-350 


(30.4±1.64)xl0"* 


15-150 


S 


>263 


>8.74 


>9.99 


6,7 


070328 


<4 


-1.11±0.04 [-2.33±0.24] 


32.6±1.24 


15-1000 


(5.9±1.2)xKr'' 


20-10000 


K/Su/S 


< 3835 


< 133 


< 690 


15 


070521 


1.35 


-0.93±0.12 


18.1±1.4 


20-1000 


(4.12±0.91)xl0"* 


20-1000 


K 


522±56 


9.22±0.71 


49.3±10.9 


18 


071020 


2.145 


-0.65±0.29 


7.71±1.36 


20-2000 


(6.04±2.08)xl0"* 


20-2000 


K 


1013±204 


8.65±1.53 


213±73 


19 


071112C" 


0.82 


-1.09±0.07 


3.0±0.4 


15-150 


8±1'' 


15-150 


S 








20 


071117 


1.331 


-1.53±0.15 


5.84±0.85 


20-1000 


(6.66±1.83)xlO-* 


20-1000 


K 


648±317 


3.59±0.52 


95.3±26.2 


21 


080319B 


0.937 


-0.86±0.01 [-3.59±0.45] 


613±13 


20-7000 


(2.26±0.21)xl0-5 


20-7000 


K 


1307±43 


142±3 


102±9.4 


22 


080319C 


1.95 


-1.20±0.10 


15.0±2.7 


20-4000 


(3.35±0.74)xl0-« 


20-4000 


K 


1752±504 


14.6±2.6 


96.1 ±21. 2 


23 


08041 3B 


1.10 


-1.23±0.25 


3.25±0.13 


15-150 


(14.0±0.58)xl0-'' 


15-150 


S 


163±34 


1.65±0.06 


14.9±0.62 


7 


080430<^ 


0.77 


-1.73±0.08 


1.17±0.05 


15-150 


(18.2±1.31)xl0-8 


15-150 


S 








7 


080602 


-1.4 


-0.96±0.63 


7.92±1.46 


20-1000 


(1.92±0.58)xl0-'' 


20-1000 


K 


>542 


>4.34 


>25.3 


24 


080603B 


2.69 


-1.23±0.64 


4.50±1.17 


20-1000 


(1.51±0.39)xl0-* 


20-1000 


K 


376±214 


9.41 ±2.45 


116±30 


25 


080605 


1.64 


-1.03±0.07 


30.2±1.2 


20-2000 


(1.60±0.33)xl0-5 


20-2000 


K 


665 ±48 


22.1±0.88 


308±64 


26 


080607 


3.036 


-1.08±0.06 


89.3±4.9 


20-4000 


(2.69±0.54)xl0"'^ 


20-4000 


K 


1691±169 


186±10 


2259±453 


27 


0806 13B 




-1.05±0.18 


22.1 ±4.1 


20-3000 


(4.76±1.31)xl0-'^ 


20-3000 


K 


733±239'' 






28 


080721 


2.591 


-0.96±0.07 [-2.42±0.29] 


88.1±7.6 


20-7000 


(2.11±0.35)xlO-5 


20-7000 


K 


1785±223 


121±10 


1038±172 


29 


080804 


2.20 


-0.72±0.04 


9.88±0.35 


8-35000 


(7.30±0.88)xlO-'' 


8-35000 


F 


810±45 


11.4±0.4 


27.0±3.3 


30 


080916A 


0.689 


-0.99±0.05 


7.02±0.21 


8-35000 


(4.87±0.27)xl0-'' 


8-35000 


F 


208±11 


0.92±0.03 


1.08±0.06 


30 


081007 


0.53 


-1.4±0.4 


1.2±0.1 


25-900 


2.2±0.2'' 


25-900 


F 


61±15 


0.17±0.015 


0.43±0.04 


31 


081121 


2.512 


-0.46±0.08 [-2.19±0.07] 


28.4±2.45 


8-35000 


(5.17±0.83)xlO-« 


8-35000 


F 


608±42 


30.5±2.6 


195±31 


30 


08 1203 A 


2.10 


-1.29±0.14 


30.5±11.2 


20-3000 


2.9±0.2'' 


15-150 


S/K 


1541 ±756 


35.0±12.8 


28.2±1.9 


32,33 


081221 


2.26 


-0.83±0.01 


27.2±0.13 


8-35000 


(24.2±0.50)xl0-'' 


8-35000 


F 


284±2 


34.6±0.2 


100±2 


30 


081222 


2.77 


-0.90±0.03 [-2.33±0.10] 


17.6±1.63 


8-35000 


(17.6±0.58)xl0-^ 


8-35000 


F 


630±31 


25.2±2.3 


94.9±3.1 


30 


090102 


1.547 


-0.97±0.01 


34.9±0.66 


8-35000 


(29.3±0.91)xl0-^ 


8-35000 


F 


1174±38 


21.4±0.4 


45.7±1.4 


30 


090201 


<4 


-0.97±0.09 [-2.80±0.52] 


67.2±5.0 


20-2000 


(7.30±1.26)xKr'' 


20-2000 


K 


<790 


<229 


<1246 


34 


090424 


0.544 


-1.02±0.01 [-3.26±0.18] 


50.1 ±1.06 


8-35000 


(9.12±0.14)xl0"'' 


8-35000 


F 


250±3.4 


3.97±0.08 


11.2±0.17 


30 



" In these cases the peak flux is in units of ph/cm^/s 

^ For GRBs without measured z the fisted peak energy refers to the observed one 
The spectrum is weU described by a power-law model and is uncostrained. 
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Table 1. Continued 



GRB 


z 


a [fi] 


Fluence 


Range 


PeakHux 


Range 


Mission 


^peak 


£■180 


^so 


Ref. 








IQ-^erg/cm^ 


keV 


erg/cm^/s 


keV 




keV 


lO'^erg 


10-^ erg/s 




090709A 


<3.5 


-0.85±0.08 [-2.7±0.24] 


91±0.7 


20-3000 


(3.9±0.6)xl0-« 


20-3000 


K 


<1341 


<229 


<442 


35 


090715B 


3.00 


-1.1±0.37 


9.3±1.3 


20-2000 


(9.0±2.5)xl0-'' 


20-2000 


K 


536±164 


21.3±3.0 


82.6±22.9 


36 


090812 


2.452 


-1.03±0.07 


26.1 ±3.4 


15-1400 


2.77±0.28° 


100-1000 


Su/K/S 


2023±663 


40.5±5.3 


96.2±9.7 


37,38 


090926B 


1.24 


-0.19±0.06 


9.81±0.16 


8-35000 


(4.73±0.28)xl0-' 


8-35000 


F 


212±4.3 


3.96±0.06 


4.28±0.25 


30 


091018 


0.971 


-1.53±0.48 


1.44±0.17 


20-1000 


(4.32±0.95)xKr^ 


20-1000 


K 


55±26 


0.80±0.09 


4.73±1.04 


39 


091020 


1.71 


-1.20±0.06 [-2.29±0.18] 


12.4±1.81 


8-35000 


(1.88±0.26)xKr'' 


8-35000 


F 


507±68 


7.96±1.16 


32.7±4.6 


30 


091127 


0.49 


-1.25±0.05 [-2.22±0.01] 


24.8±0.54 


8-35000 


(9.38±0.23)xKr'' 


8-35000 


F 


51±1.5 


1.61±0.03 


9.08±0.22 


30 


091208B 


1,063 


-1.29±0.04 


5.98±0.17 


8-35000 


(25.6±0.97)xl0-' 


8-35000 


F 


246±15 


1.97±0.06 


17.4±0.7 


30 


100615A 




-1.24±0.07 [-2.27±0.11] 


8.64±0.17 


8-1000 


8.3*0.2" 


8-1000 


F 


86*8.5* 






40 


100621A 


0.542 


-1.70±0.13 [-2.45±0.15] 


36±4 


20-2000 


(1.70±0.13)xl0-* 


20-2000 


K 


146±23 


4.35±0.48 


3.17±0.24 


41 


100728B 


2.106 


-0.90±0.07 


2.69±0.11 


8-35000 


(5.43±0.35)xl0-'' 


8-35000 


F 


404±29 


2.98±0.13 


18.7±1.2 


42 


110205A 


2.22 


-1.52±0.14 


36.6±3.5 


20-1200 


(5.1±0.7)xl0"^ 


20-1200 


K 


715±238 


55.9±5.3 


25.1±3.4 


43 


110503A 


1.613 


-0.98±0.08 [-2.7±0.3] 


26±2 


20-5000 


(1.0±0.1)xlO"5 


20-5000 


K 


572±50 


18.0±1.4 


181±18 


44 



In these cases the peak flux is in units of ph/cm^/s 
^ For GRBs without measured z the Usted peak energy refers to the observed one 
^ The spectrum is well described by a power-law model and is unconstrained. 



Table 2. Results of the statistical analysis on the £peak - £^iso and iJpeak ^ Aso correlations for (i) the total sample, (ii) the complete sample presented in this 
work and (iii) the sample of bursts belonging to the total sample but not included in the complete one (complementary sample). The table lists the Spearman's 
rank correlation coefficient (p), its associated chance probabihty, the slope and normalization of the power-law fit and the Icr scatter of the point distribution 
around the best fit line. 



Correlation 


Sample 


#GRBs 


P 


Pchance 


Slope 


Norm. 


0"sc 




Total 


136 


Q.ll 


4x10-28 


0.55±0.02 


-26.74±1.13 


0.23 


£peak ~ ^iso 


Complete 


46 


0.76 


7x10-1° 


0.61±0.04 


-29.60±2.23 


0.25 


Complementary 


90 


0.78 


3x10-" 


0.53±0.02 


-25.55±1.35 


0.25 




Total 


135 


0.74 


8x10-25 


0.49±0.03 


-22.98±1.81 


0.30 


£peak ~ ^so 


Complete 


46 


0.65 


lxlO-« 


0.53±0.06 


-25.33±3.26 


0.29 




Complementary 


89 


0.75 


3x10"" 


0.48±0.04 


-22.44±2.12 


0.30 



fined with the complete, total and complementary sample are sim- 
ilar among themselves and are about asc = 0.25 and cr^ = 0.30 
for the Amati and Yonetoku correlation. We also investigated pos- 
sible instrumental effects in the correlations found in this com- 
plete sample. We defined three sub-samples based on the instru- 
ment which derived the burst spectral parameters (see Table 1). In 
particular, we study the Amati and Yonetoku correlations within the 
sample of Swift/BAT bursts (10 events), Konus bursts (18 events) 
and Fermi/GBM bursts (13 events). The results derived from dif- 
ferent sub-samples are very similar and consistent within the er- 
rors. In fact, we foimd that the best fit slopes (normalizations) of 
the Amati correlation are 0.60±0.05 (-28.97±2.67), 0.59±0.08 (- 
28.41±4.19) and0.60±0.11 (-29.19+5.77) for the Swift, Konus and 
Fermi samples respectively. For the Yonetoku correlation we find: 
0.44±0.11 (-20.76±5.95), 0.52±0.08 (-24.61±4.49) and0.54±0.12 
(-25.87±6.37). 

Our study outlines the presence of one GRB (061021) that lies 
at 3cr limit (or more) of the tested correlations. In particular, it is an 
outlier to the Amati correlation. 

GRB 061021 has been detected by Swift/BAT (Moretti et al. 
2006) and by KonusfWIND (Golenetskii et al. 2006). It shows a 
single pulse with a duration of ~ 10 s, followed by a weak tail 



seen up to ~ 60 s. On the basis of the lag analysis (Norris & 
Barthelmy 2006) this burst is classified as a long duration event. 
The time integrated BAT spectrum is well described by a power- 
law model with photon index -1.31 ±0.06 (Palmer et al. 2006), sug- 
gesting E°ll^.> 150 keV. Golenetskii et al. (2006) report the spectral 
analysis of KonusfWIND data separately for the pulse and the tail. 
The pulse is well described by a power-law with a high-energy 
exponential cutolT, with spectral parameters a = -1.22^J| |^ and 
^peak~ ^^^-237- "^^^ ^^^^ ^^^^ modelled as a single power- 

law, with photon index -1.93+[j and a fluence which is nearly 1/3 
of the fluence of the main pulse. In order to avoid to underestimate 
the total energetics, we also include the fluence of the soft tail in 
the computation of Eiso- Nevertheless, this burst is on average ~2 
orders of magnitude less energetic as compared to bursts charac- 
terized by a similar peak energy. Note that the lack of a spectral 
analysis integrated over the whole burst duration does not affect 
the estimate of Lis;,, which is well computed by using the peak flux 
and the spectral shape of the main peak. Also in this case, the peak 
luminosity is ~2 orders of magnitude less than the average Liso of 
events with Epeak similar to GRB 061021. FinaUy, the redshift of 
tills source is z = 0.346 (Fynbo et al. 2009). 

The presence of outliers to the Amati relation is somehow ex- 
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Figure 4. Distributions of £pcak (rest frame), £iso and Ljso for bursts of 
the total (filled grey histograms), complementary (dashed dark line) and 
complete (hatched red histograms) sample (see legend in the bottom panel). 
In each panel the probability of the K-S test is also reported. 



pected from the study of the observational plane E'^^^-Fluence. 
Nava et al. (2008; 2011b) studied the distribution of BATSE and 
GBM bursts in this plane and derived respectively a 6% and 3% 
of outliers. We also stress that the definition of the outliers of the 
£^pcak - Eiso correlation refers to the present knowledge of the scat- 
ter of this correlation which is assumed to be Gaussian. It might 
well be that the scatter is still not fully known and/or that it is not 
Gaussian. The latter possibility seems to be suggested also by the 
log-normal distribution of the jet opening angles (Ghirlanda et al. 
2005) if the source of the scatter of the Epe^k - £^iso correlation were, 
indeed, the jet opening angle. 

There are 12 bursts belonging to the complete sample that 
cannot be used to directly test the correlations, due to the lack of 



measured redshift and/or £°''\ . However, we can still check their 

' peak ' 

consistency with the correlations (Fig.|2]l. We find that they are all 
consistent with both correlations. 

By taking advantage of our complete sample, for both correla- 
tions we investigate the possible evolution with the redshift of the 
best fit slopes (Fig. O. To this aim we define 4 bins of redshift, 
chosen in order to have 4 sub-samples characterized by a similar 
number of objects. We find no relation between the slope of each 
sub-sample and the redshift. All slopes are consistent with the one 
defined by the whole complete sample and this is true both for the 
£^peak - £iso and the £pcak - ^iso correlation. 

We also investigate possible differences (in terms of distribu- 
tions of the most fundamental quantities) between bursts belonging 
to the complete sample and bursts that do not satisfy the complete 
sample selection criteria. Our aim is to understand if the diff^erent 
instrumental biases present in the total sample affect in some way 
the distributions of the prompt emission properties. We find that for 
all the considered quantities (iipcak, E^o and Ljso) the distributions 
derived from the complete and from the complementary samples 
are consistent within 2a. 

Our results give support to the idea that spectral energy cor- 
relations arise because of a robust physical mechanism common 
to the vast majority of bursts. Outliers exist, but they are a few. 
Furthermore, the evidence that GRBs do evolve in redshift (Sal- 
vaterra & Chincarini 2007), while there is no sign of evolution of 
the slopes of the correlations, also suggests that the spectral peak 
energy is closely linked to the luminosity/energetics of the burst, 
another sign that selection effects are not important. 
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